Stabilization of the surface CDW order parameter by long-range 

Coulomb interaction 



P. P. Aseev and S. N. Artemenko* 
Institute for Radio- engineering and Electronics of Russian Academy of Sciences 
Mokhovaya 11-7, Moscow 125009, Russia and 
Moscow Institute for Physics and Technology 
Institutski per. 9, Dolgoprudny, Russia 
(Dated: October 18, 2011) 

Abstract 

We study theoretically formation of two-dimensional (2D) charge density wave (CDW) in a 
system of conducting chains at the surface of an insulator due to interaction of quasi ID surface 
electrons with phonons. We show that the unscreened long-range Coloumb interaction between the 
charges induced by fluctuations of the CDW phase stabilizes the finite order parameter value at 
finite temperatures, and thus the long-range order (LRO) exists. In the case of screened Coloumb 
interaction the phase fluctuations suppress the phase transition, but decay of the order parameter 
is rather slow, it obeys a power-law (A*(r)A(0)) oc r~ 7 with small exponent 7. 

PACS numbers: 71.45.Lr, 72.15.Nj, 73.20.Mf 
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In contrast to three-dimensional (3D) systems where fluctuations are usually small, in 
2D systems fluctuations can greatly affect behaviour of a system. It is well-known 1 that in 
low-dimensional (ID and 2D) continuous systems with sufficiently short-range interaction a 
long-range order (LRO) is suppressed by long-range fluctuations. These fluctuations can be 
excited with little energy cost and they are favored since they increase the entropy. Thus, 
the 2D charge-density wave (CDW) phase is also believed to be suppressed by fluctuations 
of the CDW phase. However, there are experimental evidences of existence of 2D CDW 
in layered 2 and linear-chain compounds 3 . In the latter case it was shown that the critical 
temperature for a surface CDW in NbSe^ is higher than for a bulk one. 

In a 2D superconductor which is related system, it is known 4 that there is no LRO indeed. 
Below the Berezinski-Kosterlitz-Thouless critical temperature, the correlation function of the 
order parameter obeys a power-law resulting in a pseudo-LRO in the system. 

In this paper we examine the decay of the order parameter of the incommensurate CDW 
at the surface (commensurability is expected to stabilize the LRO). 

The Mermin- Wagner theorem is applied 1 to the case of sufficiently short-range interaction. 
One of our goals is to study how the long-range Coulomb interaction may affect the LRO of 
the CDW. We consider three models of electron-electron interaction: non-screened Coloumb 
interaction, an interaction between surface electrons screened by electrons in the bulk, and 
an interaction between surface electrons screened only by surface electrons. The first case 
can be applied when the size of a sample is not greater than the screening radius. The 
second case refers to the systems where the bulk material is a semiconductor and there are 
electrons in the bulk or to systems with a gate. Finally, the third case refers to the systems 
in which the bulk material is an insulator and the only screening electrons are the electrons 
thermally excited over the Peierls gap of the 2D CDW. 

Below we set h and k-Q to unity, restoring dimensional units in final expressions when 
necessary. 

We study a system of conducting chains at the surface of a 3D insulator, electron system 
at the surface being considered as a 2D electron gas interacting with acoustic phonons in 
the bulk. We assume that a;- axis is taken along the nesting vector Q. y-axis is taken 
perpendicular to the nesting vector and parallel to the surface. Finally, z-axis is the axis 
normal to the surface. 

First, we consider the case of non-interacting electrons. The total action S is the sum 
of an action of free 2D electrons S c , action of free phonons in the bulk S p h, and action of 
electron-phonon interaction S c - P h 

S = S e + S p h + S'e-ph (1) 

Only the phonons with the wave vector component along the surface (longitudinal compo- 
nent) close to the nesting vector ±Q are relevant for formation of the CDW. If one does not 
take into account fluctuations of phonon modes one should consider only the modes with 
longitudinal component of wave vector equal to ±Q exactly and the corresponding creation- 
annihilation operators b±Q, b±Q. However, since we take into account the fluctuations we 
allow the phonons to have the longitudinal component of the wave vector that slightly differs 
from Q. We describe these phonons with components of wave-vector ±Q + q x , q y and q z 
(where q x <C Q) using creation-annihilation operators b±Q(q x ,q y ,q z ). The phonon field ip is 
introduced in a standard way 
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where u (q) = s^J(Q + q x ) 2 + q 2 + q 2 is the spectrum of acoustic phonons. Then the Mat- 
subara action of free phonons reads 



^ / uq U 2 + U 2 (q) 2 

° ph = J (2^3 Mq ^ )l ' (3) 

where T is temperature, and q stands for all the three components q x , q y , q z 

Since electrons are confined in the direction normal to the surface the electron field 
operator can be written as 

V(x,y,z) = &(x,y)w(z) 

where w(z) is the wave function of the ground mode of size quantization, w(z) decays 
exponentially in the depth of the bulk material: 

W(z) = \w(x,y,z)\ 2 = ne~ KZ (4) 

We do not take into account higher transversal modes because their contribution is expo- 
nentially small at low temperatures. Besides, the low-temperature behaviour of the system 
and particularly the properties of CDW are determined by the electrons with energy close 
to the Fermi energy, and, therefore, with momentum close to ±Q/2. Thus we can represent 
the electron field as 

^(a:, y) = e iQx/2 Jj + {x, y) + e- iQx/2 ^.{x, y), 

where ip±(x, y) vary smoothly in comparison with the correspondent exponential factor. The 
Matsubara action of free electrons reads 

S c = J ^- 2 dr {r + [d T - e + ] ^ + p_ [d T - e_] V- } , (5) 

where e± = Eo(k ± Q/2) and Eo(k) is a spectrum of free electrons. 
The term describing the electron-phonon interaction reads 

= - J rfVrf 3 rrfr^;(r)^_(r)x 

xW(z)g(r-r')ip Q (r') + c.c. (6) 

We assume that the Fourier transform of interaction potential g(k) = J d 3 rg(r)e~ tkr may 
depend on the momentum k. Although we do not specify the particular dependence, note 
that such a dependence appears if one considers an electron-phonon interaction in a system 
with anisotropy. 

It is convenient to introduce the order parameter as follows: 

A(fc x , k v ) = J ^¥ Q , qz ( k x, k y )g qz (k x , k y )W qz (7) 

where W qz is the Fourier transform of the square of transversal wavefunction (4) W qz = 
J dzW(z)e~ igzZ . The electron-phonon contribution to action (6) can be rewritten in terms 
of A in a quite simple form 

S c - ph = - J d 2 rdr {^-A + A*CV+} ( 8 ) 
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If one does not take into account fluctuations of the phonon modes then A does not 
depend on coordinates and can be taken real. In this case one can calculate electron Green 
functions in imaginary time 



g±± = <v4v>±> = 

and obtain the following expressions 

G + 4e,k) = G_ + (e,k) 
G ++ (e,k) = -G*__ = - 



Se—Sc — ph 



e 2 + e + A 2 



e 2 + e + A 2 



where £ = £+( - fc - > 2 £ — ~ Vpk x . These Green functions correspond to an electron spectrum 
with a gap e = a/£ 2 + | A| 2 . This result is similar to 3D case 5 . 

It is convenient to integrate out the phonon fields and thus to derive an effective action 
for A. Using the definition for A (7) we obtain 



J Vip* Q V<p Q e- s »» = J Vip* Q V<p Q VA*VAV/3*Vf3x 

x exp {-S ph + (3*(k x , k y ) [A(k x , k y )- 



<P*Q,q,( k x> k v)9* qz (kx, k y ) w qf^ 



+ c.c. 



where P(k x ,k y ), /3*(k x ,k y ) are Lagrange multipliers for A* and A. Performing Gaussian 
integration over (p,(p*,/3,/3* we obtain an effective action for A 



S A = -TJ2 



dkxdky \A(u,k)\ 2 /F(u,k) 



F(u,k)= J 



(2tt)= 

dq z ujj(k,q z ) 
2ir u 2 + Uo(k, q z 



\9l(h,k y )\\wl 



(9) 
(10) 



We can expand F(uu, k) at small u and k 

F(^k) = 9 -t 



u 2 + s' 2 k 2 



s 2 (Q 2 - K 2 ) 

where go = g(k = 0), and s' can be roughly estimated as 



12 _ Jg(Q)-9o\ 



9o 



(11) 



We consider the case when Q > k. If Q < k the formation of CDW due to interaction with 
the phonons in the bulk is impossible. 
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Now we can minimize the total action given by (1), (5), (8), (9) and thus find the classical 
solution for the order parameter A (i.e. the solution that does not take into account fluc- 
tuations of phonon field and, therefore, fluctuations of the order parameter). The equation 
for the classical A resembles the self-consistent condition for the gap in the BCS theory 



A 



/ 



= tanh —-de, 
2 2T 



where A = 4 ^°* a , a y is a lattice constant in the direction normal to the nesting vector. 
The cut-off parameter at high energies is of order of Bp and thus the solution at T = is 
A = Epe~ x . The result is similar to the case of CDW in 3D. The main feature of 2D systems 
is that in constrast to 3D, fluctuations are important in 2D. The main contribution to the 
fluctuations of the order parameter is given by fluctuations of the phase x where A = A e %x . 
The amplitude of the order parameter determines the Peierls energy gap, and fluctuations 
of the ampltitude are described by the mode with a gap in the spectrum, in contrast to 
the long wavelength fluctuations of the phase which are related to gapless modes, so that 
at q — > they do not affect the total energy. Thus at T < A we can ignore amplitude 
fluctuations, and it is convenient to rewrite the total action in terms of x Then the phonon 
action reads 

d 2 q X(u,q)x(-u,-q) 



S A = -T\A fJ2 



(2k)' 



F 2 (u,q) 



The total effective action is 



S CS [X*,X] = S A -lnJ Vr + Vi> + Vr-Vi>-e 



—S e —S ( 



e~ de— ph 



(12) 



Expanding the exponent in powers of x an d leaving the second order terms, we perform 
Gaussian integration over ip±, ip±. After some algebra we obtain the following effective 
action 



s cS = tJ2 



d 2 k 2IAI 2 



■x 



(2vr) 2 g 2 K 

A u 2 + s' 2 \k\ 2 \u 2 + v 2 k 2 x 



2 s 2 (Q 2 + k 2 ) 6 |A| 2 

+ 



X(e,k)x(-e,-k) 



A A u 2 + v\k 2 x 



12 |Anl 2 



Green function for x can be easily calculated 

(x(u,k)x(-uj,-k)) = 
(x(u,k)x(u,k)) = 



V X X 2 (oo,k)e' s ^ 
n\v 2 J(Sv F ) 



uj 2 + v Ik 2 + s' 2 k 2 



(13) 
(14) 



where v 2 ~ v'p ' V J A | 2 ' V ; is a velocity of excitations along the x-axis. Note that it follows 
from (11) that the velocity of excitations s' along the y-axis is non-zero if we take into ac- 
count a dispersion of electron-phonon interaction g{k). However, one would obtain non-zero 
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velocity of excitations along the y-axis if one took into account non-ideal nesting conditions. 
Thus we can assume that in a real setup the value of s' is of order of sound velocity s. The 
correlation function for phase x can De calculated as 



(x(r,t)x(r + 6r,t)) = 



= E 



' H{ (x(w, k)x(-w, -k)) e ikSr — 



(2vr) 2 

vrA f <fk e* fr ygg + gg 

32 VF 7 (2 7 r) 2 v ApJ+^ 2T 1 bj 

If the temperature T = then the correlation function for phase is finite and therefore 
the LRO exists in the system. Otherwise, if the temperature T > 0, the integral in (15) 
diverges at small k when 5r = 0. Thus there is no long-range order at finite tempera- 
tures. However, the correlation function for the order parameter (A*(r, t)A(r + 5r,t) = 

| Ao |2 e -([ X (5r-,t)- X (0,t)]x(0,t)) obeyg the power _l aw 

(A*(0)A(r)) = |A | 2 (rQr-, 7o = ^^ 

The exponent 7 is of order of ^, where 6 is the Debye temperature. 

In the presence of Coulomb interaction there is an addition contribution to the total 
energy, and thus to the total action (1), due to an interaction of charges induced by phase 
fluctuations. Given the phase Xi the electron density can be calculated as 



e 



-s 



Thus we obtain the electron density induced by phase fluctuations 5p = —ed x x/ . 
The corresponding term in the action responsible for Coulomb interaction is given by 



Sc = J dTd 2 r 1 d 2 r 2 Sp(r 1 )V(r 1 - r 2 )5p(r 2 ) 



We consider three different models for Coulomb potential V{r\ —r 2 ): non-screened Coloumb 
potential, Coloumb potential screened by free electrons in 3D material, and Coloumb po- 
tential screened by thermally excited electrons in 2D electron layer. 

First we consider non-screened Coulomb potential V(r) = e 2 /(e*r), where e* is effective 
dielectric constant. If the surface is a boundary between two media with dielectric constants 
ei and e 2 , then e* = (e x + e 2 )/2 Calculating the Green function for x by using (13) we obtain 

\nvl/(8v F ) 

(x(w,q)x(-w,-Q))- 



2 + 2^ii g + s >2 2 
e*a y vp 



The singularity at low frequencies is now integrable if T > as well, so the mean square of 
fluctuations of phase (x{ r ^)x{ r ^)) is finite and therefore the LRO exists in the system. 

Consider Coulomb interaction screened by free electrons in 3D, for example by electrons in 
a gate or by conduction electrons in a semiconductor. The interaction is described by Yukawa 
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potential V(r) = r / rD , where td is a screening radius. The Fourier transformation of 
this potential reads 



V(q) = j d 2 rV(r)e- iqr = 2ire 2 / [ey q 2 + r^ 2 ) , 
and the Green function for x calculated using (13) has the form 

Xnv 2 J(Sv F ) 



uj 2 + v 2 x [l + l2\^)q 2 x + s' 2 q 



y 



The integral over u and q diverges at T > as in the case when there is no inter-electronic 
interaction, thus there is no LRO. However, the correlation function for A obeys a power-law 
with different exponent. 



<A*(0)A(r)> = A c 



-73 



73 = 7o 



12e*\hv F 




The exponent 73 7o, and although there is only a short-range order in the system, the 
correlations decay more slowly and the system behaves like it has a real LRO provided the 
screening radius is large enough. 

Finally, we consider the Coulomb interaction screened only by electrons thermally excited 
over the Peierls gap in the conducting 2D layer of width a z . The Fourier transformation of 
the potential is given by 

V{q) = 2TTe 2 /[e*{\q\ + a z /rl)] 

Density of screening electrons and, therefore, the screening radius r D depends exponentially 
on temperature oc n~ l oc e A//T . 

The Green function for x calculated using (13) reads 



(x(u,q)x(-u,-q)) 



\7iv 2 J(8v F ) 



cu 2 + v 2 x ( 1 + 12A-g--fM q 2 x + s' 2 q 2 

x V e*VF a y a z 1 lx 1 y 

The integral over u and q diverges at T > again, and the correlation function for A obeys 
a power-law but with a smaller exponent 



(A*(0)A(r)> = A (ra z /r 2 D ) 



-72 



72 = 



16^3 s' Hv F 



D 



The exponent 72 73 70 for similar values of screening radius, and the correlations 
decay at greater distance than in the case of 3D screening or non-interacting electrons. 
Strictly saying, there is no LRO in the system but due to the very slow power-law decay of 
correlations there is a pseudo-LRO. 
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